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Abstract 

Aero-assisted orbit transfer vehicles have the potential for significantly reducing the fuel 
requirements in certain classes of orbit transfer operations. Development of a nonlinear feedback 
guidance law for performing aero-assisted maneuvers that accomplish simultaneous change of all 
the orbital elements with least vehicle acceleration magnitude is discussed. The analysis is based on 
a sixth-order nonlinear point-mass vehicle model with lift, bank angle, thrust and drag modulation 
as the control variables. The guidance law uses detailed vehicle aerodynamic and the atmosphere 
models in the feedback loop. Higher-order gravitational harmonics, planetary atmosphere rotation 
and ambient winds are included in the formulation. Due to modest computational requirements, 
the guidance law is implementable on-board an orbit transfer vehicle. The guidance law 
performance is illustrated for three sets of boundary conditions. 

Introduction 

The concept of using maneuvers in the planetary atmosphere for orbital plane change was first 
advanced by London [1], This concept did not receive significant attention in the literature until the 
early eighties. Since then there has been significant research activity in evaluating the payoffs and 
in synthesizing guidance laws for the aero-assisted orbit transfer mission. For instance. References 
2 and 3 examines the performance payoffs in employing such maneuvers. Guidance law 
development using various vehicle model approximations are discussed in References 4-10. A 
Lyapunov optimal feedback control law for the orbit transfer mission was proposed in Reference 
11. Using an approximate guidance law. Reference 12 discusses the effect of earth rotation on the 
aeroassisted orbit transfer mission. 

The aeroassisted orbit transfer maneuver guidance laws reported in the literature [2-12] achieve a 
desired velocity vector orientation at the point of exit from the atmosphere while optimizing a 
specified performance index. The optimization objectives have included the minimization of energy 
loss during the atmospheric pass, minimum fuel, minimum aerodynamic heating, and minimum 
flight time. Minimum acceleration magnitude criterion is used as the performance index in the 
present research. This criterion will ensure that the emerging guidance law results in the least 
structural load on the vehicle, permitting implementation on-board the orbit transfer vehicle. 

The present paper addresses the problem of simultaneously controlling all the six components of 
the vehicle state vector at the atmospheric exit. Such a formulation enables simultaneous in-plane 
and out-of-plane aeroassisted orbit transfer maneuvers together with a control of the flight time, 
latitude and longitude at the exit from the atmosphere. Control of flight time is useful for 
establishing the time-of-perigee passage or the argument of perigee of the resulting orbit, while the 
control of terminal latitude and longitude are important to assure adequate ground station visibility. 



The focus of the present research is on the development of a feedback guidance law for on-board 
implementation. Feedback linearization theory will be used to develop an aero-assisted orbit 
transfer guidance law that achieves a specified terminal state vector in a specified time period. This 
paper is an outgrowth of the preliminary research outlined in Reference 13. 

The feedback linearization concept has been previously applied to various nonlinear guidance 
problems [14, 15] with significant success. The present research differs from the studies given in 
References 4 - 12 in several ways. For instance, a planet centered Cartesian coordinate system is 
employed in the present formulation. The effects of planet oblateness, atmosphere rotation and 
ambient winds are included. Note that neglecting the atmosphere rotation will result in an 
overestimate of aerodynamic forces in prograde maneuvers, and will lead to underestimating the 
aerodynamic forces in retrograde orbit transfer maneuvers. The present analysis allows for bank 
angle, lift, drag modulation and thrust control. While the use of lift, bank angle and thrust control 
have been previously considered, the drag modulation is proposed here to aid in satisfying the 
flight time constraint The formulation uses a sixth-order nonlinear vehicle model with four control 
variables. Assuming that all the state variables can be measured or accurately estimated, and that 
the system nonlinearities are known with reasonable accuracy, this model is transformed into a 
linear, time-invariant form. The guidance problem is then formulated and solved in the feedback 
linearized coordinates. The resulting solution is subsequently transformed back to the original 
coordinates to obtain the nonlinear feedback law. The nonlinear guidance law incorporates all 
known details about the vehicle and ambient conditions. For instance, it can utilize real-time 
estimates of atmospheric density and winds if these were available. Consequently the present 
guidance law can be expected to have an acceptable performance even while operating under off- 
nominal conditions. Further details on feedback linearization, guidance law development, and 
numerical evaluation are furnished in the ensuing sections. 

Point-Mass Spac ecraft Model 

The present analysis will only consider the atmospheric portion of the aeroassisted orbit transfer 
mission. Thus, the following development assumes that the thrust impulse required for entry into 
the planetary atmosphere has been applied. Point-mass equations of motion for a vehicle flying 
over an oblate spheroidal planet in a planet centered inertial frame is given by the vector differential 
equation 
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x, y, z are the vehicle position components measured in the planet centered inertial frame (X, Y, Z) 
shown in Figure 1. The X-axis of the inertial frame points along a chosen principal direction, 
while the Z-axis is aligned with the planet's rotation axis. The Y-axis completes the right handed 
triad. In the case of earth, the principal direction is often chosen as the direction of the Vernal 
Equinox [16]. Figure 1 also shows the topocentric inertial frame (Xj, Yy, Zt) used for the 
calculation of the flight path angle y and the heading angle %. In equation (1), the quantity T is the 
vehicle thrust, D is the vehicle drag, L the lift, m the vehicle mass and <J> the bank angle. 0 is the 
longitude with respect to the principal direction and \| / is the latitude. Vehicle mass m and thrust T 
can be functions of time. g x , gy, gz are the components of acceleration due to gravity. In all that 

follows, a dot over variables will indicate differetiation with respect to time, and the superscript T 
will denote the matrix transpose operation. 



Fig. 1. Coordinate Systems 


Given the vehicle position in the geocentric inertial frame, the corresponding latitude and longitude 
can be calculated as: 
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( 2 ) 

Longitude with respect to any specified references on the planet such as the Greenwich meridian in 
case of earth can be found by subtracting the instantaneous longitude of the reference meridian with 
respect to the Vernal Equinox direction from the longitude 0. If the vehicle velocity vector with 
respect to the geocentric inertial frame is known, the flight path angle y and the heading angle x 
can be computed as follows. First compute the velocity vector in the topocentric system using the 
transformation: 
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The Xt axis of the topocentric frame points in the southern direction, the Yj axis points in the east 
direction and the Zj axis points in the direction of local vertical, see Figure 1 for details. The flight 
path angle with respect to the local horizontal and the heading angle with respect to south can then 
be computed as: 


= tan T-l, y tan 1 , W - 
luj |_Vu +v z _ 


(4) 


If desired, the heading angle with respect to true north can be computed as (re — %)• Figure 2 shows 
the topocentric frame together with the vehicle flight path axis system. 



Fig. 2. Topocentric Frame and the Flight Path Axis System 
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Figure 2 also shows the sign convention employed for the definition of the vehicle bank angle. The 
vehicle thrust and drag acts along the Zf axis, while the vertical component of lift acts along the - 
Xf axis. The horizontal component of lift acts along the -Yp axis. It is assumed here that the thrust 
acts along the velocity vector and that the vehicle employs coordinated maneuvers. In cases where 
these two assumptions are not admissible, an additional transformation incorporating the angle of 
attack and the angle of sideslip needs to be included in the analysis. 


In order to compute the aerodynamic forces acting on the vehicle, the vehicle relative velocity with 
respect to the atmosphere needs to be found. The atmospheric rotational velocity with respect to the 
inertial frame and the local wind velocities can be included to yield the components of the vehicle 


velocity vector 


relative to the atmosphere as 
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where the radial position of the vehicle with respect to the planet center is given by 


r = ^/x 2 + y 2 + z 2 

( 6 ) 

0) is the planet's rotational rate and u w , v w , w w are the ambient wind velocity components defined 
in the topocentric frame. These wind components can be specified as a function of latitude, 
longitude and altitude. The vehicle altitude at any time instant can be computed as 


h=^ 


x 2 + y 2 + z 2 — r c 


(7) 

The variable r c is the distance between the planet surface and the origin of the coordinate system 
measured along the direction of r. The vehicle velocity with respect to the planet atmosphere is 
given by 

V.=V?+^+wf 


( 8 ) 

This velocity should be used in the calculation of aerodynamic forces. The gravitational 
acceleration model used in the present analysis can include any number of harmonics. For instance, 
the gravitational model may be of the form [16]: 
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Here \i is the gravitational constant for the planet, r is the radius vector, r e the planet radius, and 
J2, J 3 ... are the harmonics of the planet's gravitational potential. 

The vehicle aerodynamic lift and drag are given by the expressions 

L = ipV 3 sC L , D = ipV. ! s(C D 

( 12 ) 

The drag coefficient Cd could be a function of the lift coefficient, Mach number and altitude. The 
control variable CDm is the component of drag that can be modulated using aerodynamic devices. 
Cl is the lift coefficient, another control variable in the problem. The bank angle <J> and the vehicle 
thrust T are the additional control variables in the problem. The air density p can be a function of 
altitude, latitude and longitude. The variable s is the reference area used in the calculation of 
aerodynamic forces. 

The vehicle model elaborated in the foregoing is nonlinear and time varying. In the present form it 
is useful only for numerical studies. In the next section it will be shown that the model can be made 
amenable to analysis through a series of transformations. 


The spacecraft nonlinear point-mass model may be transformed to a more convenient form by 
assuming that all the state variables required for the calculation of the right hand side of equation 
(1) are available from measurements. In this case, the right-hand-side can be treated as a three 
component pseudo-control vector, yielding the equations of motion as: 
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Ui, U 2 , U 3 are the pseudo-control variables. These quantities are unknown at this stage because 
although the state variables and the system nonlinearities are known, the control variables are yet to 
be determined. Physically, these quantities are the spacecraft acceleration components in the 
geocentric inertial frame. If the pseudo control components Ui, U 2 > U3 wer c available, the actual 
control variables T, Cl, CDm, <t> can be computed using an inverse transformation as follows. 
First compute the vehicle flight path angle, heading angle, latitude and longitude using the 
measured vehicle velocity and position components. The force vector [Fx, Fy, Fj"^" in the flight 
path axis system is next computed using the expression: 
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The vehicle lift coefficient can then be calculated as : 


2^ + F, 2 

PV. 2 s 




The bank angle <j> is given by 


0 = tan '[F y /F 2 ] 


(15) 


(16) 


Next, the thrust and drag modulation coefficient can be computed based on the sign and magnitude 
of the force component F z as follows: 
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Thus, if the pseudo-control variables Ui, U 2 , U 3 are known, the real control variables T, Cl, 
CDm, and <t> can be recovered using the inverse transformation (14) - (18). 
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In the case of pure aerobraking maneuvers as in Reference 1, 2, 4 - 12, vehicle thrust and drag 
modulation capabilities may not available. In this case, the inverse transformation has to account 
for this fact by considering the third component of the vector equation (14) as a constraint of the 
form : 
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This constraint can be u»ed to determine one of the pseudo control variables in terms of the other 
two. The vehicle lift coefficient and bank angle can then be calculated as indicated. However, the 
guidance objectives have to be modified in this case because the simultaneous control of flight 
time, position and velocity control at the atmospheric exit will no longer be possible. Such a 
formulation will not be examined in the present research. In the ensuing, it will be assumed that all 
the four control variables are available. 

Guidance L aw Development 

The previous section dealt with an approach for making the nonlinear vehicle model amenable to 
analysis. In this section, the guidance problem will be formulated in the transformed coordinates 
and solved. Inverse transformation of this solution to the original coordinates will then yield the 
nonlinear feedback guidance law. The following development will assume that all the state 
variables required for computing the feedback law are known or can be estimated. Once the 
guidance problem is solved with perfect information, the effects of incomplete or imperfect 
information can be investigated using this solution. 

It is assumed here that altitude, flight path angle, heading angle, latitude, longitude at atmospheric 
exit and the time of flight through the atmosphere are all specified. Since the present formulation 
uses a planet centered inertial frame, these conditions can first be translated to the Cartesian frame 

using well known relationships in astrodynamics [16]. For instance, given the flight path angle 
Yf, heading angle Xf> latitude yf, longitude 0f and the velocity at atmospheric exit Vf, the 

components of the terminal velocity vector in the planet centered inertial frame can be calculated 
using the expression 
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( 20 ) 


The objective of the present research is to synthesize an implementable guidance law that achieves a 
specified set of final conditions. As mentioned elsewhere in this paper, the performance index is 
assumed to the integral of the square of the vehicle acceleration. Primary objective is to synthesize 
a maneuver gui dan ce law that results in least structural load on the vehicle. Since the magnitude of 
a vector is invariant with respect to rotational transformations, the vehicle acceleration magnitude is 
the same as the magnitude of the pseudo control variables. In addition to being mathematically 
convenient, this performance index may also lead to lower energy loss during the maneuver. With 
this, the aeroassisted maneuver guidance problem may be defined as: determine the vehicle lift, 
bank angle, drag modulation coefficient and thrust that achieve 

min \ f[U* + U’+U*ldt 
u„u„u, I J L 1 


( 21 ) 


given the position vector and velocity vectors at the atmospheric entry and exit points, subject to 
the differential constraint (13). The final time tf is assumed to be specified. If desired, flight time 
can be included in the performance index as in Reference 15. 


In order to solve this optimal guidance problem, the variational Hamiltonian [17] is first defined as: 


H = — [Uf + + U*]+ A,x + A 2 y + A 3 z + X 4 \J l + A 5 U 2 + A 6 U 3 
2 


( 22 ) 

(23) 
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The initial and terminal conditions on the costates are unknown. Expressions (24) can be integrated 
together with (23) to yield 

A 4 — Cj ■“ A,t, Aj c 2 A 2 t, X 6 c 3 A 3 t 

(26) 

ci, c 2 , C 3 are constants of integration. Substituting in (25), the control variables become 


The costate equations and the optimality conditions for this problem can be obtained as 

Aj = A 2 — A 3 = 0 
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(27) 


Substituting (27) in the feedback linearized equations of motion and integrating results in: 
x = x 0 + t|^A 1 ^-c 1 j, y = y 0 + t^ 2 |-c 2 j, z = z 0 + 1|\ j - c 3 j 

(28) 

X = x 0 + x 0 t + y^ 1 |-c 1 j, y = yo+yot + j[^2"C 2 J z = z 0 + z 0 t + 

(29) 

The constants cl, c2, c3, Xi, X2, X3 are the unknowns in these equations. These may be calculated 
from the given terminal conditions. Using the specified exit velocity components and the final time, 
the costates Xi, X2, X3 can be computed as 


Xj — — [x f x 0 +c,t f ], X 2 - [y f y 0 +c 2 t f ], X 3 - [z f z 0 +c 3 t f ] 

l r l f 

(30) 

These may next be substituted in the equations (29), together with the specified terminal 
conditions yield 
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At this stage, an open-loop maneuver guidance law can be obtained by substituting for ci, c2, C3, 
Xl, X2, X3 in (27). It can be verified that the guidance law satisfies the strengthened Legendre- 
Clebsch and Jacobi necessary conditions [17]. As a result, all the trajectories emerging from the 
present analysis afford a strong local minimum. Further, since there exists a unique optimal 
trajectory for every specified set of boundary conditions, these trajectories provide a global 
minimum [18]. 


The guidance law can next be converted to feedback form by assuming that the current time is the 
initial time. In this case, the difference between the current time and final time is the time-to-go 
parameter, denoted by the variable tg 0 . The resulting guidance law in terms of pseudo-control 
variables may be substituted in the inverse transformation to yield the nonlinear guidance law for 
aeroassisted orbit transfer as: 
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where 



(33) 

The vehicle lift coefficient, bank angle, drag modulation coefficient and thrust can be calculated 
from the left hand side of this equation as in (15)-(18). Note that all the quantities on the right hand 
side can be computed from the current and the specified final vehicle states, and time-to-go. The 
performance of this guidance law will be illustrated using a point-mass simulation of a aero- 
assisted orbit transfer vehicle in the next section. 

Results and D iscussions 

The guidance law developed in the previous section is next implemented on a point-mass model of 
a hypothetical maneuverable re-entry vehicle shown in Figure 3. The performance data for this 
vehicle [6] is given in Table 1. As presently configured, this vehicle has no drag modulation 
capabilities. However, for the purposes of guidance law evaluation it will be assumed, that such a 
capability exists. Additionally, it is assumed that the orbit transfer maneuvers occur with earth as 
the primary. 


11 




Fig. 3. Maneuverable Re-entry Research Vehicle 


The vehicle simulation featured a fourth order Runge-Kutta integration scheme with a choice 
between linear or exponential interpolation of the U. S. standard atmosphere up to a maximum 
altitude of 300000 feet Four gravitational harmonics were included in the simulation. 

Vehicle Mass 335.67 Slugs 

Reference Area 125.84 ft^ 

Zero-Lift Drag 0.032 

Coefficient 

Induced Drag 1 .4 

Coefficient 

Maximum Lift 0.151 

Coefficient 

Table 1. Maneuverable Re-entry Vehicle Data 

Three sets of boundary conditions are considered in this paper. The First one is that of a pure 
orbital plane change maneuver, the second one a simultaneous change in both in-plane and out-of- 
plane parameters, and the last one being an in-plane maneuver. Right duration for all the three 
cases were assumed to be 400 seconds. The trajectory boundary conditions used in the study are 
summarized in Table 2. 


12 



CASE.l 



Latitude, 

Longitude. 

Altitude, 

Flight Path 

Heading 

Velocity. 


deg. 

deg. 

feet 

Angle. 

Angle. 

feet/s 





deg. 

deg. 


INITIAL 

0 

0 

200000 

-0.148 

90 

25946 

FINAL 

10 

20 

200000 

0.0 

130 

14000 




CASE 

JL 




Latitude. 

Longitude. 

Altitude, 

Flight Path 

Heading 

Velocity, 


deg. 

deg. 

feet 

Angle, 

Angle. 

feet/s 





deg. 

deg. 


INITIAL 

0 

0 

200000 

-0.148 

90 

25946 

FINAL 

10 

20 

200000 

10 

130 

14000 




CASE 3 




Latitude. 

Longitude. 

Altitude. 

Flight Path 

Heading 

Velocity. 


deg. 

deg. 

feet 

Angle, 

Angle, 

feet/s 





deg. 

deg. 


INITIAL 

0 

0 

200000 

-0.148 

90 

25946 

FINAL 

0 

20 

200000 

10 

90 

14000 


Table 2. Trajectory Boundary Conditions 

In the first orbit transfer maneuver scenario, the vehicle enters the atmosphere over the equator 
heading due east The vehicle travels north and exits the atmosphere at 10 degree north latitude and 
20 degrees north longitude. The desired exit flight path angle is zero. The only difference between 
the first and second case is the exit flight path angle. A rather high terminal flight path angle in 
Case 2 corresponds to a highly eccentric target orbit. In Case 3, the exit flight path angle is same is 
Case 2. However, the terminal heading angle and latitude are different. 

During the course of these simulations it was found that the linear interpolation of the atmosphere 
produced comers in the control history. For this reason, all the results given in this paper have 
been generated using an exponential atmosphere model. Figure 4 illustrates the vehicle altitude 
history for the first set of boundary conditions. The vehicle initially descends to an altitude below 
175000 feet before climbing back to meet the terminal altitude specification. The vehicle latitude- 
longitude evolution is given in Figure 5. The guidance law exactly meets the specified terminal 
altitude, longitude and latitude. 
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210000 



0 100 200 300 400 

TIME, 3 

Fig. 4. Altitude Vs Time (Case 1) 



The temporal evolution of the vehicle inertial velocity is shown in Figure 6. The corresponding 
flight path angle and heading angle are shown in Figures 7 and 8. Note that the heading angle with 
respect to North is shown in Figure 8. The heading angle history is not monotonic because the 
vehicle is required to simultaneously meet all the components of the terminal velocity and position 
vectors. 
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Inertial Velocity, ftls 



Fig. 6. Vehicle Inertial Velocity Vs Time (Case 1) 
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Time seconds 

Fig. 8. Heading Angle Vs Time (Case 1) 

(Measured With Respect to North) 

The lift coefficient corresponding to this trajectory is shown in Figure 9. From this figure it can be 
observed that the vehicle tends to use a higher lift coefficient at higher altitudes. This strategy is 
consistent with the specified objective of minimizing the vehicle acceleration magnitude since larger 
lift coefficients at lower altitudes tends to not only increase lift, but also contributes to drag. 

0.12 

„ 0.10 
A 

A 

£ 0.08 

a 

4) 

S 0.06 
S 

J 0.04 
0.02 

0 100 200 300 400 

Time, seconds 

Fig. 9. Lift Coefficient Vs Time (Case 1) 

The vehicle bank angle history is shown in Figure 10. As the vehicle begins the atmospheric 
maneuver, it banks more than 90 degrees to the left and decreases the bank angle to zero before 
beginning a bank to the right. If desired, the positive bank angle component can be eliminated by 
an appropriate re-definition of the terminal latitude and longitude. 
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100 



Fig. 10. Bank Angie Vs Time (Case 1) 

The vehicle drag modulation coefficient along the trajectory is illustrated in Figure 1 1. The vehicle 
does not call for any drag increase during descent, but requires it as the vehicle is climbing out to 
meet the specified terminal velocity and position conditions. 



Fig. 11. Drag Modulation Coefficient Vs Time (Case 1) 

The corresponding thrust history is illustrated in Figure 12. The thrust is used for velocity control 
during the vehicle descent. It is important to recall here that the present study achieves the control 
of all the six orbital parameters. As a result, the thrust and drag modulation histories given in 
Figure 1 1 and 12 depend on the specified boundary conditions. 
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Fig. 12. Vehicle Thrust Vs Time (Case 1) 

The vehicle Mach number history along the orbit transfer maneuver trajectory is illustrated in 
Figure 13. The vehicle starts the maneuver at about Mach 26 and has a nearly 14 Mach at 
atmospheric exit. 



Fig. 13. Mach Number Vs Time (Case 1) 

In Figures 14 through 21, compares the trajectories emerging from the three sets of boundary 
conditions given in Table 2. The altitude histories for the three cases are compared in Figure 14. 
Lowest altitude corresponds to the case of simultaneous change in all the orbital elements. Due to 
the high terminal flight path angle requirement, Case 3 also requires significant descent into the 
atmosphere. 
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300000 

| 200000 
% 

1 

g 100000 

0 

0 100 200 300 400 

Time, seconds 

Fig. 14. A Comparison of Altitude Histories 

The vehicle inertial velocities are compared in Figure 15. The velocity histories are nearly identical 
for case 1 and 2, while the maneuver corresponding to case 3 evolves at a much lower average 
velocity. The Mach number histories given in Figure 16 reveal that although the inertial velocities 
for case 1 and 2 were nearly identical, the vehicle experiences different airspeeds and densities due 
to the differences in the altitude histories. 

26000 
" 24000 

s 

£ 22000 
| 20000 
£ 18000 
| 16000 
3 14000 

12000 

0 100 200 300 400 

Tim* seconds 

Fig. 15. A Comparison of Inertial Velocities 
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30 



0 100 200 300 400 

Time, seconds 

Fig. 16. A Comparison Mach Number Histories 

The flight path angle along the three maneuvers are shown in Figure 17. The vehicle lift coefficient 
histories given in Figure 18 reveal that a much lower lift magnitude is used for in-plane changes 
when compared with out-of-plane changes. 



Fig. 17. A Comparison of Vehicle 
Flight Path Angles 
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Fig. 18. A Comparison of Lift Coefficient Histories 

The corresponding bank angle histories are given in Figure 19. Of the three histories given in 
Figure 19, the one corresponding to case 3 is the most interesting one. Upon atmospheric entry, 
the vehicle initially assumes an inverted position which lasts for about 40 seconds. Subsequently it 
uses zero bank angle during the remaining duration of the maneuver. 



0 100 200 300 400 


Time, seconds 

Fig. 19. A Comparison of Vehicle Bank Angle Histories 

Finally, the vehicle drag modulation and thrust histories are given in Figures 20 and 21. Both case 
2 and 3 call for significantly higher thrust than Case 1. These increased thrust magnitudes are used 
primarily for offsetting the increased drag caused bu flight at lower altitudes. On the other hand, 
both these cases use lower drag modulation than Case 1. 


21 





Fig. 20. Comparison of Drag Modulation Coefficients 



Fig. 21. A Comparison of Thrust Histories 

Conclusions 

This paper presented the development of a nonlinear feedback guidance law for the aero-assisted 
orbit transfer maneuvers. It is useful for the simultaneous change of both in-plane and out-of-plane 
orbital elements. The analysis employed a nonlinear point-mass model of the spacecraft. Other 
modeling details were the inclusion of arbitrary-order gravitational model, atmosphere rotation, 
ambient winds and an atmosphere density model of arbitrary complexity. 
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For guidance law development, the vehicle model was first transformed into a linear, time invariant 
form. In order to minimize the structural load on the vehicle, the performance objective was 
assumed to be the minimization of the square of the vehicle acceleration. The velocity vector 
orientation and the position vector components at the atmosphere exit were specified. The guidance 
problem was solved in feedback form in the transformed coordinates. Inverse transformation of 
this solution produced the nonlinear guidance law for aero-assisted orbit transfer. The guidance 
law is in closed-loop, state feedback form and uses the vehicle performance data in the feedback 
loop. Numerical results using a maneuverable re-entry vehicle data were given. 
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Appendix - A 

Implementation of the AOTV Guidance Law on a Macintosh Computer : 

All the results given in this report were generated using a point-mass simulation implemented on 
a Macintosh computer. In this section, some of the details of the Graphical User Interface for the 
AOTV simulation are given. 

Four menu items at the top of these figures are used to control the program operation. The three 
choices available under the TRAJECTORY menu items are : Data, Run, Terminate. 
Selecting the Data item results in the screen display: 

TRAJECTORY CONSTRAINTS SUBSYSTEMS OUTPUT 


ID SIMULATION DATA, SECONDS „ ^ - 

STEP SIZE, S 

| 1.000 


FINAL TIME, S 

100.000 


OUTPUT INTERUAL, S 

2.000 





These edit fields can be used to change integration step size, final time and the output interval. 
The other two choices in this menu causes the simulation run or termination. 

The CONSTRAINTS menu items allow for the input of boundary conditions, and time- 
acceleration weights in the performance index. As an illustrative example, the following 
illustrates the screen format for the "boundary conditions' menu item choice. 
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TRRJECTORY CONSTRAINTS SUBSYSTEMS OUTPUT 


TRRJECTORY BOUNDARY CONDITIONS 


INITIAL COND. 
FINAL COND. 


INITIAL COND. 
FINAL COND. 


ALTITUDE 


LATITUDE LONGITUDE 


200000.000 

.000 

.000 

200000.000 

20.000 

20.000 

UELOCITY 

FLIGHT PATH 

HERDING 

259456.703 

-.148 

.000 

1400.000 

.148 

40.000 


OK 


The SUBSYSTEMS menu item allows the input of model subsystem data such as: 
Aerodynamics, Atmosphere Density Model, Wind Profile, and Gravitation Harmonics. The 
screen displays for Aerodynamics, Atmosphere Density Model, Gravitation, and 
Wind Profile menu items are shown in the following figures. .. 
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r 


TRAJECTORY CONSTRAINTS SUBSYSTEMS OUTPUT 


=r~i 3 mb UEHICLE AERODYNAMIC DATA I 



UEHICLE MASS, Slugs 

| 335.674 


REFERENCE AREA, ft*ft 

125.840 


ZERO-LIFT DRAG COEFFICIENT.... 

.032 


INDUCED DRAG COEFFICIENT 

1.400 


MODULATION DRAG COEFFICIENT 

.001 


MAHIMUM LIFT COEFFICIENT... 

.151 



m 

k 


Clicking on the OK edit field will cause the program to accept the values typed in the edit fields. 
As currently configured, the user interface accepts data up to four decimal places. 
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r TRAJECTORY CONSTRAINTS SUBSYSTEMS OUTPUT 

jgliMsi atmospheric densitv models 


O EHPONENTIflL (i) PIECEUJISE LINEAR 



At the present time, the program permits the use of either an exponential or piecewise linear 
model of Earth. Future versions may include the atmospheric data for various other planets also. 

S imilar ly, the program currendy incorporates a multi-harmonic gravitational model of earth. The 
number of gravitational harmonics included in the formulation can be increased or decreased by 
clicking on 'up' or 'down' fields. Inclusion of the gravitational acceleration of other planets will 
be a future item. 
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r TRAJECTORY CONSTRAINTS SUBSYSTEMS OUTPUT 


ID_jlB number of grauitationrl harmonics i 


NUMBER OF HARMONICS UP) 

4 

EM 






WM 


The ambient wind model currently accepts the wind velocity and its azimuth specified as a 
function of altitude. Inclusion of latitude-longitude dependence will be undertaken in the future. 
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Finally, the OUTPUT menu item allows access to the simulation results, either as an output 
data file compatible with plotting package such as CRICKET™ or as interactive data display. 

At the time of this writing, the program kernel and the graphic user interfaces have been tested in 
a stand alone mode. Integration and comprehensive testing of these packages will be of future 
interest. 
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